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We present recent results for the prediction of the Bg—B® lifetime difference from lattice Heavy Quark Effective 
Theory simulations. In order to get a next-to-leading order result we have calculated the matching between QCD 
and HQET and the two-loop anomalous dimensions in the HQET for all the AB=2 operators, in particular for 
the operators which enter the width difference. We present results from quenched and, for the first time, from 
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unquenched simulations. We obtain for the B s —B g lifetime difference, Fb 3 = (5.1 ± 1.9 ± 1.7) x 10 and 
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— — = (4.3 ± 2.0 ± 1.9) x 10 from the quenched and unquenched simulations respectively. 



1. Introduction. 

In the Standard Model, the width difference 
( ) in the B®— B® system is expected to be the 
largest among bottom hadrons and could be mea- 
sured in the near future. The Standard Model 
prediction for r ^ a relies on an operator product 
expansion, where the short distance scale is the 
b quark mass The theoretical expression for 
the width difference reads, schematically [0, 
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K(G(z) + G s (z)TZ(m b ) + 8 1/m ) 



(1) 



where if is a factor which encloses known pa- 
rameters, z — ml /ml, G(z) and Gs(z) are NLO 
Wilson coefficients calculated in Jjl and 5\/ mb are 
0(l/mb) corrections which depend on matrix el- 
ements of higher dimension operators (see |l| for 
details). All the non-perturbative contribution, 
at leading order in l/mj, comes from 1Z, 



Tl(m b ) 

where 

O l = 
O s = 
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~ (B s \0 L (m b )\B s ) : 
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We present results for the 1Z parameter and for 
other AB=2 B-parameters from lattice HQET in 
the quenched approximation and also, for the first 
time, from unquenched simulations. We have per- 
formed the matching between lattice HQET and 
continuum QCD at NLO. Until now this match- 
ing was known at NLO only for the Ol operator 
(see ref. 0). We have extended this calculation 
to all AS = 2 operators by computing the NLO 
matching between QCD and HQET in the con- 
tinuum and the relevant two-loop anomalous di- 
mensions ||. 

2. B-parameters from the lattice. 

In the HQET there are four independent opera- 
tors which we write in terms of the B-parameters, 
defined as vacuum insertion deviations: 



B(p), (4) 



(B d{s )\0 L (p)\B d{s )) 

l f B d{s ) M k (s ) 
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= B s (p) (5) 



(B d{s )\0 LR (p)\B d{s )) f 1+ 2 A T 1 



■)LR 



<M»»( 1 + 6 ^,-- s B- W (T) 



•i P M 

3JB d ,y U B d(s ) 



K =0.1580 K =0.1570 L_=7 t,/a=3 



where, 



X = 



Ml 
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(nib + m d ( s) ) 2 
Ol and Os are denned in eq. (||) and 
LR = 1^ (1 - 75 )g6 7M (l+ 75) q, (9) 

0£ fl = 6(1-75)96(1+75)?. (10) 

The lattice B-parameters are extracted from 
the large time behavior of the ratio between three- 
and two-point correlation functions (see H): 



Cp^-hM) ti,fa-joo (B a |Oi(a)|B a ) 
C(-*i)C(t 2 ) ~ ' \(0\A (a)\B s )\ 2 



where 



C Q% (h,t 2 )=J2 {0\Ao (x! , hWO, 0)4 (f 2 , ta) |0) 
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and is the HQET axial current. 

The results presented here are obtained from 
two simulations: a quenched one performed by 
APE coll. with the SW-Clover action on a 24 3 x40 
lattice with 600 gauge configurations at f3 = 6.0 
with aT 1 = 2 GeV and an unquenched one carried 
out by T%L coll. (gauge configurations) and APE 
(correlation functions) with the Wilson action on 
the same volume with 60 gauge configurations at 
P — 5.6 with two degenerate sea quarks corre- 
sponding to a- 1 = 2.51, 2.54 GeV. See refs. § 
and (7) respectively for details. 

In figure 1 we show the plateaux for the four 
ratios Roi in the unquenched case. We observe 
good plateaux over large time-distances for all op- 
erators (In ref. the corresponding figures in 
the quenched case were presented). 

In [|| we pointed out that different groups using 
different methods found in quenched simulations 
that the lattice B-parameters in the static limit 
are very close to 1. As one can see from figure 1, 
this result holds also in unquenched simulations. 

3. 1/rrib dependence of the B-parameters. 

Having performed the calculation described in 
the previous section one obtains the B-parameters 
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Figure 1. The ratios i?o 4 for different operators 
computed at k va i = 0.1570 and k sea = 0.1580 for 
ti/a = 3 and with a smearing size of 7. 



in QCD to leading order in 1/mj,. By using vac- 
uum insertion approximation (VIA) for the sub- 
leading operators || and the fact that the bare 
lattice B-parameters are very close to 1, we have 
shown |l(J that all AB = 2 B-parameters, de- 
fined as vacuum insertion deviations, have small 
0(A/m,b) corrections (A = Mb — rrib), 



Bi{m b ) = Bi(m b ) 



O( 0.3 A) +0(-L), (11) 
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where Bi(mb) and Bi(m,b) are the full and the 
static QCD B-parameters, respectively. The fac- 
tor 0.3 is an estimation of the deviation from VIA 
in the calculation of the subleading operators and 
0(a s ) corrections. 

4. Results. 

Before presenting our results, let us stress that 
m s and nib parameters in previous sections are 
the corresponding pole masses. The latter coin- 
cides with the expansion parameter of the HQET 
because we have set the residual mass to zero. We 
calculate the b quark pole mass from the running 
MS mass, which can be accurately determined. 
Since our computation is performed at NLO, we 
use the perturbative relation between the pole 
and the running mass at the same order. From 



Table 1 



Values of the B-parameters in 
scheme (see text). 


the NDR-MS 


B-Parameters 


Clover 


Wilson 


quenched 


unquenched 


B(m b ) 
Bs(m b ) 
B LR (m b ) 
Bl R (m b ) 


0.83(5)(4)(5) 
0.96(8)(5)(5) 
0.94(5)(4)(5) 
1.03(3)(5)(5) 


0.74(3)(8)(5) 
0.87(2) (11) (5) 
0.98(4)(6)(5) 
1.10(2)(9)(5) 



Table 2 

Values of B, B$ and 1Z from different groups. 



Grc 


urn 


B(m b ) 


Bs(mb) 




[1 




0.85(3)(11) 


0.80(2)(10) 


-0.91(5)(17) 


[2]J 


0.91(3)1°;°° 


0-86(2)±8:g 


-0.93(3)1°;°° 


This work 
quenched 


0.83(5)(6) 


0.81 (7)(6) 


-0.95(7)(9) 


This work 
unquenched 


0.74(3) (9) 


0.74(2)(10) 


-0.97(5) (15) 



the world average running mass [ fLl| , m b {m b ) = 
4.23 ± 0.07 GeV one obtains m b = 4.6 ± 0.1 GeV 
Notice that the contribution of m s is very small 
because it always appears divided by m b . 

In table [l] we present the values for all AB=2 
B-parameters. The first error in table 1 comes 
from lattice simulations and includes statistical 
and systematic errors. The second one is an es- 
timate of the error due to the uncertainties in 
the values of the lattice coupling constant and 
to higher order contributions to the matching. 
The third is an estimate of l/m b corrections 
to the static result (see eq.(|il"|)). The analy- 
sis of the perturbative matching is explained in 
detail in our previous work where we gave 
B(m b ) ( - quenched ) = 0.81 ± 0.05 ± 0.04. The tiny 
difference is due to the fact that there we used, in 
the perturbative evolution, a number of flavours 
rif = 4 instead of rif = as in the present paper. 
From table 1 we see that the systematic errors 
are larger in the unquenched case than in the 
quenched case. This is due to the fact that the 
renormalization constants for the Wilson action 
are bigger than those for the SW-Clover one, and 
consequently, the contribution of higher orders in 
a s are more important. Taking into account the 
errors in table 1 we cannot establish any clear 
unquenched deviation from the quenched result. 

The B-parameters in table [j] are QCD scheme 
dependent, they all are in the NDR-MS scheme. 
To make contact with the computation of ||, we 
have subtracted the evanescent operators in the 
renormalization of Ol and Os as in pL For B LR 
and Bg R we use the prescription of |{L3| . 

For 1Z (see eq.(2)), which contains all the non- 
perturbative contribution to the width difference, 



at leading order in l/m b , we obtain 

TZ(m b )^ ue ^ = -0.95(7)(9) (12) 

TZ(m b ^ un ^ = -0.97(5)(15) (13) 

where as usual, the first error comes from lattice 
simulations, and the second is systematic due to 
the uncertainty in the perturbative matching and 
to the contribution of higher orders in l/m b . 

5. Comparison with other recent results 

In table |^ we present our result for the B- 
parameters relevant for r Bs compared with 
other recent quenched determinations. In [0] the 
Bs is defined as in eq. (@) but in terms of the run- 
ning mass instead of the pole mass. Therefore, 
we have multiplied our value by (m b (m b ) / m b ) 2 
in order to compare with their result. On the 
other hand, in the definition used in JTJ] does 
not appear the factor X (see eq. (||)), therefore 
we have divided their result by X (m b /m b (m b )) 2 . 
Notice that in spite of using different methods to 
obtain the B-parameters, there is a good agree- 
ment between the three quenched computations, 
in particular in the value of the ratio TZ. In the 
unquenched case the central values of B and Bs 
are lower but still compatible within errors. Nev- 
ertheless, the ratio TZ, is in very good agreement 
with the quenched computations. 

From equation (Q) we get our prediction, 
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= (4.3 ± 2.0 ± 1.9) X 10" 



(15) 



The first error is systematic obtained from the 
spread of values of all input parameters in eq. (0) . 
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Figure 2. Linear and quadratics fits of Bg(mi>) in 
1 /Mb to the ref. data (open circles) combined 
with the static point (filled square). The filled 
diamonds correspond to physical B s meson mass. 



The second one comes from the uncertainty in 
the value of 5i/ mb . Since in the estimate of Si/ mb 
the operator matrix elements were computed us- 
ing VIA and the radiative corrections were not 
included, we assume an error of 30% [§. As can 
be seen, this parameter is still affected by a large 
uncertainty, so that a precise determination of the 
width difference requires the computation of the 
subleading matrix elements using lattice QCD. 
This simulation is missing to date. 

Our result is to be compared with the present 
experimental status |15fl 
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(16) 



As can be seen the central values are rather dif- 
ferent but still compatible within the large errors. 

Finally, in figures 2 and 3 we present different 
fits in 1/Mb to the quenched data of refs. PJl^| 
for Bs and the renormalization group invariant 
parameters Bs d and Bb s - As we pointed out in 
sec. 3, the HQET predicts that the 1/Mb cor- 
rection to the static B-parameters is small and 
the contribution of the quadratic term cannot be 
neglected. To study this dependence, we perforin 
different fits: a linear fit, a quadratic fit, including 
the linear term, and a pure quadratic fit including 
and excluding the static point. The conclusion is 
that the quadratic fits give a smaller value for the 
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Figure 3. The same of figure 2 for Bb s and Bb s 
with the data of ref. p6l. 



physical B-parameters than the linear one and are 
in better agreement with the static value. Nev- 
ertheless, the predictions of all quadratic fits are 
still compatible within errors with the linear fit. 
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